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Abstract 

The purpose of the paper is to find explicit formulas describing the joint distributions of the first 
hitting time and place for half-spaces of codimension one for a diffusion in M""*"^, composed of one- 
dimensional Bessel process and independent n-dimensional Brownian motion. The most important 
argument is carried out for the two-dimensional situation. We show that this amounts to computation 
of distributions of various integral functionals with respect to a two-dimensional process with inde- 
pendent Bessel components. As a result, we provide a formula for the Poisson kernel of a half-space or 
of a strip for the operator (/ — A)"/^, < a < 2. In the case of a half-space, this result was recently 
found, by different methods, in pj. As an application of our method we also compute various formulas 
for first hitting places for the isotropic stable Levy process. 

1 Introduction 

Bessel processes appear in many important theoretical as well as practical apphcations. A systematic 
study, in the frame of one-dimensional diffusions, was initiated by H.P. McKean in [TB]. Bessel processes, 
in particular, play an important role in a deeper understanding of the local time of the standard Brownian 
motion (see Ray-Knight Theorem in e.g. [20] )• Their generahzations serve as various models in applica- 
tions ^Oj. In our presentation, when considering Bessel processes, we follow the exposition of [20], Ch. 
XI. 

The purpose of the paper is to exploit a correspondence between harmonic measure of the operator 

in 

(/_A)°/2, 0<a<2, 

for a half-space or a strip and the joint distribution of the first hitting time and place for appropri- 
ate sets of codimension one for the (n -|- l)-dimensional Bessel-Brownian diffusion, that is, the process 
Y{t) = where is one-dimensional BES^ "/^^ process and i?" is independent 

n-dimensional Brownian motion. Although we do not appeal to Molchanov-Ostrovski representation (see 
[19]) our results exemplify their principle that some jump Markov processes (here: isotropic stable and 
relativistic Levy processes) can be regarded as traces of the appropriate Bessel-Brownian diffusions. With 
this regard we also mention the paper P] where the eigenvalue problem for the Cauchy process in i?" 
was transformed into a kind of mixed eigenvalue problem for the Laplace operator in and also the 

work [7] which proposes to study n-dimensional non-local operators by means of (n-|- l)-dimensional local 
operators. 

The paper is organized as follows. In Prehminaries (Section 2) we enclose some essentials about 
various special functions, indispensable in the sequel. We also provide a basic information about squared 
Bessel processes, Bessel processes and their hyperboHc and trigonometric variants, needed in the sequel. 
A detailed discussion of generalized Bessel processes is carried out in Appendix. 
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In Section 3 we state our basic observation that the m-harmonic measure for some regular sets £) C 
for the n-dimensional Q-stable relativistic process can be read from the joint distribution of the first exit 
time and place for the set D C M"+^ (with the trace D) for our Bessel-Brownian diffusion. The proof of 
this, intuitively apparent result, is postponed to Appendix. When m = we obtain the same statement 
for the usual harmonic measure for the standard isotropic a-stable process. 

In Section 4 we first compute the joint distribution of the first hitting time and place for the vertical 
positive axis in the two-dimensional space for the process Y. Next, we deal with the same process Y 
but this time hitting two horizontal half-lines (— oo, —1] and [1, oo). This is crucial for further purposes. 
We rely on stochastic calculus; in particular, we apply appropriate random change of time and compute 
various integral functional of Bessel processes and generalized Bessel processes. The resulting formulas 
are, in the first case, in terms of modified Bessel functions; in the second one - in terms of spheroidal 
wave functions. In the end of this section we apply our results to provide a purely probabilistic method 
of computing the Poisson kernel of the interval [—1, 1] for the standard isotropic a-stable process. 

In Section 5 we generalize previous results to multidimensional case. In the view of our principle 
from Section 3, the first part of Section 5 yields the formula for the Poisson kernel of a half-space for the 
operator (/ — A)°/^, < a < 2, recently found in |Hj, using rather analytical methods. The second part 
gives a description of the Poisson kernel for a strip; also for the standard isotropic a-stable process. 



2 Preliminaries 

In this section we collect some preliminary material. For more information on modified Bessel functions, 
Whittaker's functions and hyper geometric functions we refer to PP and [H]. For questions regarding 
Bessel processes, stochastic differential equations and one-dimensional diffusions we refer to PD] and to 

2.1 Special functions 
Modified Bessel Functions 

Various potential-theoretic objects in the theory of the relativistic process are expressed in terms of 
modified Bessel functions and K^. For convenience we collect here basic information about these 
functions. 

The modified Bessel function Lq of the first kind is defined by (see, e.g. [TT], 7.2.2 (12)): 

^-M = (|)'E(|)" \,r(t^, + i) . (1) 

where € M. The modified Bessel function of the third kind is defined by (see [TT], 7.2.2 (13) and (36)): 
K^z) = [/_^(z) - /^(z)] , ^iZ, (2) 



n e Z. (3) 



We will also use the following integral representations of the function K^{z) (|llj, 7.11 (23) or pjj, 8.432 
(6)): 

[■CO 2 

K^{z) = 2-^-^z^ e-'e-^t~^^-Ut, (4) 
Jo 

where 3?(z^) > 0, | argz] < |. 

In the sequel we will use the asymptotic behavior of as a function of real variable r: 

w = wj^^r^Y^ ^-o+, (5) 



r(?? + i) V2 
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where g{r) = f{r) means that the ratio of g and / tends to 1. When r ^ 1 we have 

I^(r) = (27rr)i/2ell + i?(r)], (6) 

where E{r) = ©(r^-*^), r — > oo. 

Confluent hypergeometric function and Whittaker's functions 

The Laplace transforms of some additive functionals of Bessel processes are given in terms of the Whit- 
taker's functions. We introduce some basic notation and properties of these functions. 
The confluent hypergeometric function is defined by 

$(a,6;z) = 5]ff 1^, 6/0,-1,-2,..., 

fc=0 ^"^'^ 

where z is a complex variable, a and h are parameters. Here (a)^ = r(a + /c)\r(a) denotes the Pocchamer 
symbol. For | arg < vr and 6/0, ±1, ±2, ... we define a new function 

^' ^) = vJ^lr^K^ ^) + ^-HtV^ ^(1 + o - 6, 2 - 6; z) 

i[\ + a — 0) i (aj 

called the confluent hypergeometric function of the second kind. The Whittaker's function of the first and 
the second kind are defined respectively by 

M,,^(z) = 2^^+1/2 ^(1/2 -Av + /x,2/i + l;z), (7) 
W^^^{z) = z''+^/^e^'/'^^{l/2- K + fi,2fi + l;z), |argz|<7r. (8) 

The confluent hypergeometric function of the second kind satisfies the relation 

^{a,b;z) = z^^^^{l + a-b,2-b;z), |arg2;|<7r, 

which implies that 

W^,^{z) = W^,.^{z). (9) 
The following asymptotic formula ([I], 13.5.8 p. 508), valid for 1 < 5R(6) < 2, 

^ia,b;z) = ^^z'-'' + Oil), \z\^0 

together with the definition of the Whittaker's function give that for < /i < 1/2 we have 

= ^-"^'^-^'( r(v''!l,) ^"- + 0(i)). M^o. (10) 

Hypergeometric function and Legendre functions 

For c / 0, —1, —2, ... we define the hypergeometric function by 

77/1, ^ (a)n(6)n „ 11^1 

2Fi{a,b;c;z) = }_^ ^, z , \z\<l. 

n=0 



If K(c - a - 6) > then 



^/ , ^ r(c — a — 6)r(c) 

2F1 a,6;c;l = -) ' \ . 11 

1 (c — ajl (c — oj 
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The function 2^1(0, b; c; z) is a solution to the hypergeometric equation 

z{l - z)y"{z) + [c-{a + b+ l)z]y'{z) - ahy{z) = 0. (12) 

If c is not an integer then the other independent solution is given by 

z^-^ 2^1(0 + 1 - c, 6 + 1 - c; 2 - c; z) = z^~^{l - zY'^'^ 2-^1(1 - a,l - b-2 - c; z). (13) 

The Legendre functions are solutions of Legendre 's differential equation 

(1 - z^)y"{z) - 2zy\z) + \v{v + 1) - - z-'Y^Uz) = 0. (14) 

Making appropriate substitutions it can be reduced to the hypergeometric equations and consequently 
the Legendre junction of the first kind is defined by 

and the Legendre function of the second kind is 

Qiiz) = e-2-^- VV^I(^±ii±ll.-^-.-Hz^ - 1)-/^ 2F,(^ + 1, ^^±i^; u + 

2 

The Wronskian of the pair {Plf{z),Q',!^{z)} is equal to 



T{u+l) ' ' ' 2 ' ' 2 ' ' 2'z2' 



(16) 



1 ^....2. r(^ + i)r(^) 
i-z^ r(i + ^)r(i±p^)' 



(17) 



2,2 Bessel processes 

Time change of squared Bessel processes 

The basic material concerning Bessel processes is taken from [20], Ch. XI. 
We begin with a definition of a square of Bessel process: 

Definition 2.1. Let P{t) be a Brownian motion. For every 6^0 and x ^ 0, the unique strong solution 
of the equation 

dZ{t) = 2y^\Z{t)\dp{t) +ddt, Z{0)=x, 

is called the square of 6- dimensional Bessel process started at x and is denoted by BESQ^{x). 6 is the 
dimension of BESQ^. The square root of BESQ^{a^), a ^ 0, is called the Bessel process of dimension 5 
started at a and is denoted by BES^{a). 

We introduce also the index v = (5/2) — 1 of the corresponding process, and write BESQ^'^^ instead 
of BESQ^ if we want to use instead of 5. The same convention applies to BES^'^\ 
For V > —1, the semi-group BESQ^'^\x) has a transition density function 

qt\x,y) = \(^Y'\-^-+yy^'l,{^/t) ior x>0 (18) 

and 

^Y\^, y) = (2t)-'^-it-('^+i)r(i/ + l)-i y'^e-^/^t ^ ^^^^ 

ly denotes here the modified Bessel function of the first kind. It is well known that for —1 < 1/ < the 
point is reflecting; for ^ the point is polar. 

The infinitesimal generator of BESQ^ is equal on C^(0, 00) to the operator 
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If X is a BESQ\x), then for any c > 0, the process c^^X{ct) is a BESQ\x/c). 

Next, we begin with introducing the two dimensional process Y = (11,12) with independent compo- 
nents as Yi being the squared BESQ^~"^'^^ process and Y2 - the standard Brownian motion. Equivalently, 
Y = (11,12) satisfies the following stochastic diferential equations 



clYi = 2y^\dBi + (2 - a)dt 
dY2 = dB2, 



(21) 



where {Bi,B2) denote the standard two dimensional Brownian motion {EBi{t) = 0, EBf{t) = t for 
i = 1, 2), and where < a < 2, li(0) = yi ^ and 1^2(0) = 2/2, with 2/2 S M. It is well-known that Yi ^ 0. 
Consequently, the absolute value in the first equation of ( [2T| ) can be discarded. Also, for — 1 < < 0, the 
point is (instantaneously) reflecting and the process li is (pointwisely) recurrent. 

Let {Bi,B2) be the standard Brownian motion in independent from {Bi,B2). We consider another 
pair {Xi,X2) of independent squared Bessel processes BESQ^~°'/'^'> defined by the following system of 
stochastic differential equations 



dXi = 2y/X~idBi + {2 -a) dt 
dX2 = 2^/X2dB2 + {2-0) dt, 



(22) 



such that Xi(0) = xi, ^2(0) = X2, where xi,X2 ^ 0. 

We define a function / = (/i, /2) : ^ by /(xi, X2) 
Using Ito formula we obtain 



(4xij;2,a;2-a;i)- Let {Zi, Z2) = f{Xi,X2). 



dZ, 



dZ2 = 

Observe that if we put 



4X2^X1 + 4XidX2 

4X2 2^fx[dBi + 4X2(2 - a)dt + 4Xi 2^/X2dB2 + 4Xi(2 - a)dt 
4^X1X2(2^2 dBi + 2y/XidB2) + 4(Xi + X2)(2 - a)dt 
2v^(2v^d5i + 2^/x[dB2) + 4{Xi + X2){2 - a)dt, 
-2Jx'idBi+2JX2dB2. 



dUi 
dU2 



2^/X2dBl + 2^/Xl dB2, 
- 2 Jx~idBi + 2 JX2dB2 , 



then the process (f/i, U2) is a continuous martingale starting from such that 



d{Ui){t) 

d{U2){t) 

d{Ui,U2){t) 
If we define the following integral functional 



Mt) 



4{Xi + X2)dt, 
4(Xi + X2)dt, 
0. 



iXiis) + X2is))ds 



and its generalized inverse function 



cji(s) = inf{t > : Ai{t) > s}, 



then we get (see [E] Theorem 7.3 p. 86) that the process {Ui{ai{s)) , U2{(Ji{s))) is a standard 2- 
dimensional Brownian motion (/3i,/32)- Consequently we obtain 

dZioai = 2y/\Zi o ai\dPi + {2 - a)dt, 
dZ2 oai = dP2- 



This means that Y = Zoai, where Y = (11,^2) is the process defined by (|2l[). 
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Time change of generalized Bessel processes 

In this section we consider the process X = (Xi, X2) given by the following stochastic differential equations 



dXi 
dX2 



,J\l-Xl\dBi-'^-^Xidt 
2 - a, 



^\Xl-l\dB2 + ^^1X2! dt, 



(23) 



such that Xi{{)) = xi, X2(0) = X2, where \xi\ ^ 1 and 3:2 ^ 1. Here -Bi, B2 are two independent standard 
Brownian motions in M. It follows that that for |Xi(0)| ^ 1 we have ^ 1 for all t > 0. The 

boundary points 1 and (—1) are instantaneously reflecting. 

Analogously, X2{t) ^ 1, whenever X2(0) ^ 1 and the boundary point 1 is also instantaneously 
reflecting. For justification of these statements, see Appendix. 

The first process is a version of Legendre process; the second one - of hyperbolic Bessel process. 

Obviously the processes Xi, X2 are independent. Moreover, the generators of the processes are given 

by 



G2 



l-x^ d^ 



2- a d 

~ir''di' 



X 



2 dx^ 
2 - 1 2- a d 

2~d^ ^ ~^''d^' 



(24) 
(25) 



respectively. 

Let h = (/ii,/i2) : — > by hi{xi,X2) = (1 — ^D^x^ — 1) and h2{xi,X2) = xiX2- We define the 
process Z = {Zi, Z2) = h{Xi,X2). Using Ito formula we get 



dZi 



-2 Xi {Xl - 1) dXi + 2 X2 (1 - Xl) dX2 - (X| - 1) (1 - Xl) dt + (X| - 1) (1 - Xf) dt 



= 2^/\Z^\{-Xi ^Xl -ldBi+ X2 - Xl dB2) + (2 - a)(X| - Xj) dt, 
dZ2 = X2dXi + XidX2 



= {X2^l-XldBi + Xi^Xl-ldB2) 
Observe that for 

dWi = 

dW2 = X2^l - Xl dBi + Xi^Xl - 1 dB2, 
we have d{Wi){t) = d{W2){t) = X| - Xl and d{Wi, W2){t) = 0. Thus for 



Xl JX^ -ldBi + X2 \ 1- Xl dB2 



{Xl{s)-Xl{s))ds 



(26) 



and (T2 = inf{t > : A2{t) > s} we get that {Wi{a2{t)),W2{a2{t))) is a standard two dimensional 
Brownian motion (/?i,/32) and consequently we have 

dZi o a2 = 2\/ Z\ o 02 dPi + (2 — a)dt, 
dZ2 o (72 = dP2- 



Comparing this with (21 ) gives Y = Z o a2. 
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3 Relativistic stable processes and Bessel-Brownian diffusions 

Assume that < a < 2. A Levy process X'^{t), living in M", is called the a-stable relativistic stable 
process with parameter m > if its characteristic function is given by 

^Og(^,x'-(t)> ^ gmtg-i(|z|2W/-)°/2 ^ zeR"". 

For m = we obtain the standard isotropic a-stable process. The generator of X"^{t) is given by ml — 
(m2/"I- A)"/2. For more information about relativistic processes we refer the reader to [5] and [22] . 

For an open set l) C M" we define = inf{t > 0; X™'{t) ^ D} . The X-harmonic measure of the set 
D is defined as 

where A is a Borel subset of M". In this paper we are interested only in the case A = m and we will denote 
P^'™ as P^. Note that the relativistic process killed at an independent exponential time with expectation 

1/m has the generator equal to — (m^/"/ — A)°/^. Therefore can be regarded as the harmonic measure 

of D for the operator —{m?/°^I — A)"/^. If this harmonic measure is absolutely continuous with respect 
to the Lebesgue measure on we call the corresponding density the m-Poisson kernel of D and denote 
byP^(5,-). 

Next, let Y{t) = {Yi{t) B"^ {t)) be an (n + l)-dimensional diffusion with independent components, 
where Yi(t) is BES'^-"'/^'> and B'^it) = {B2{t), B^it), . . . , Bn+i{t)) is the standard Brownian motion in 

The following proposition exhibits the fact that finding m-harmonic measures is equivalent to finding 
some hitting distributions of the process Y{t). 

Proposition 3.1. Let D C M" be open. Assume that F = W^\D has the interior cone property at every 
point. Let D = M"+i \ ({0} x F). Let x = (0,S) G {0} x l) . Define td = inf{t > 0; Y{t) ^ D} and 
assume that P^{t£) < oo). Then for every Borel A C M" we have 

P^{x,A) = E^[e-'^^°;B''{Tn) e A]. 

The conclusion is valid also for m = 0, that is for the harmonic measure for the standard isotropic a-stable 
process. 

The proof of the above proposition is provided in the Appendix. 

4 Hitting distributions in 

We begin this section with considering two dimensional case, which is crucial for further computations in 
higher dimensions. 

4.1 Hitting distribution of a positive vertical half-line in 

We compute here the joint distribution of the first hitting time and place of the positive vertical axis 
for the process {r\~°'^'^\ B{t)) with independent components, where r[~"'^'^'^ is a BES'^~°-/'^'> process and 
B{t) is the standard Brownian motion. We always assume here that < a < 2. Define the following two 
sets 

H = {(xi,X2) G [0, cxd) X [0, oo) : xi > 0}, 
Di = {(yi,y2) eMx [0,oo) : (yi = 0) A(y2 >0)r. 

Let denote by td^ the first exit time of the process {r[~°'^'^\ B{t)) from the set Di, i.e. td^ = inf{i > 
: (i?S""/^\p(t)) ^ Di}. Observe that the BES^~'^/'^\a) process pj""/^^ hits the exactly when 
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BESQ^-'^/^\a^) process {R 



i-a/2).2 



hits the 0. Therefore, we consider, as in (21 1 processes (11,12) with 



Yi = (ii(-°/2))2 Y2 = B. We thus have 

TD^ = inf{t > : Y{t) ^ Di}. 

We denote by tr the first exit time of the process X = {Xi,X2) with independent BESQ^~°'^'^^ compo- 
nents, as defined by (22), from the set H, i.e. 

TH = inf{t > : X{t) i H] = mi{t > : Xi{t) = 0}. 



Recall that f{xi,X2) = {4x1X2, X2 — xi). 




Figure 1: Simulated paths of X = {Xi,X2) and image Y = (li,l2) under mapping /. 

Lemma 4.1. The distribution of {t£)-^,Y{ti:)^)) with respect to P^y^^^y^) is the same as the distribution of 
{Ai{th), f{X{TH))) with respect to {xi,X2), where f{xi,X2) = (yi,y2)- 

Proof. We recall that Ai{t) = 4 /q(Xi(s) + X2{s)) ds. It is easy to see that f{H) = Di and that the 
function / is bijective on H. Moreover, from the time change property proved in the previous subsection 
we get 

f{X{a^m=Y{t), 

where the equality in distribution is meant for the underlying processes. Hence in order to prove the 
lemma we may and do assume that f{X{ai{t))) = Y{t), t > 0. We have 



TDi = inf{t > 
= inf{t > 
= inf{t > 



Y{t) i Di} 

i H} 

= inf{y4i(n) > : X{u) ^ H} 
= Ai{th) 

and 

Y{td,) = f{X{a^{TD,))) = /(X(ai(Ai(rj,)))) = /(X(th)) = {0,X2{th)). 

□ 

The proof of the main result of this subsection is based on well-known formulas describing distributions 
of some integral functionals of quadratic Bessel processes, stated here in the following two lemmas. For 
convenience of the reader we include proofs. 
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Lemma 4.2. Let X be a BESQ^'^^ process with u > —1. Then for x > the following holds: 



exp (-y X{s)ds^ ■,X{t) G dr 



For x = we obtain 
exp 



A 



X{s)ds ■,X{t) G dr 



2sinh(tA) 



1 



^ coth(tA) T 



xrX 



sinh(tA) 



1/ 



r^X 



r(l - u) (2sinh(tA))'^+i 



g-^coth(a) 



Proof The above follows directly from the formula for the distribution of the corresponding integral 
functional for the quadratic Bessel Bridge BESQ'f\x,y) which can be found, e.g. in [2^3, Ch. XI, 
Corollary 3.3, p. 465: 



exp 



A2 

' 2 



X{s)ds 



A 



sinh(A) 



exp 



^^(l-Acoth(A))'^' 



sm"h(A) ) 



where Q^,y^ is the distribution of the quadratic Bessel Bridge BESQ^^\x,y). 

The above formula, the scaling properties of BESQ'^'^^ processes and the formula (18) for the density 
of BESQ^"^ give the first formula of the Lemma. The second formula is obtained from the first one by 
limiting procedure as 6 ^ 0, taking into account ([S]) □ 

Lemma 4.3. Let X be a BESQ^"^ process with u > —1. Define ti, = inf{s : X{s) = b}. Then for 
X > b > we have 



u.y{x) = E^[exp{—jTh 
For —l<u<0 and b = we obtain for x > 

A2 



£;'^[exp(-7ro 



TO 



X{s)ds)] 



p/ i^+l+A/2 N 



x{>^+i)/2r(\u\y 



(27) 



Proof. By the form of the generator of the process X (|20j) and the general theory of Feynman-Kac 
semigroups [9] we infer that the function u^{x) satisfies for x > 6 the following diferential equation 



2 + 2{v + 1) - {{\^x)/2 + 7) u,{x) = . 



dx'^ 



dx 



(28) 



The two linearly independent solutions of (28) are of the form 



^(x) = x-(''+i)/2m_^/2a,|h/2(Ax) and </)(x) = x-^'+^^'^W_^/2x,\,\,2{\x) , 



where M_ 



7/2A>|/2 



and VF_^/2A,|!/|/2 ^tre Whittaker's functions defined in ([7|) and ([Sj). Taking into account 



-//^^ii"!/^ -r/^^in/^ O l_f 

the fact that the function u^{x) is bounded and u^{b) = 1 we obtain the first formula. The second formula 



follows from the first one and the asymptotic formula (10) which gives 

r(|H) 



6-(-+i)/2T^_,/2A,|H/2(Afc) = A(i-^)/25-e-^V2 
Hence, we obtain 



r( 



i^+l+A/2 



-A'^6'^ + 0(l) 



lim b-^'+^^/^W. 



■7/2A,|i/|/2 



(A6) 



r(IH) 



which proves the second formula. 



□ 
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Theorem 4.4. For (ijf, 5(0)) = {zi, Z2) £ Di, zi > and r > we have 



_£;(2i,z2) 



e~^^«i;-B(rDj G dr 



{\z\ + Z2)^(k| - Z2)"^ 
2Tr(f)W4 



For zi = and Z2 = u < we get 



sm(7ra/2) /-n\"/^ e'^^^"") 



TT 



r — M 



(29) 



Proof. Let (/) G C^(0,oo) and A > 0. From Lemma [41] we obtain that 



e-^^«i0(y(ri5j) 



It is convenient to carry out our computations for A/2 instead of A. From the independence of the processes 
Xi and X2 and the fact that th is determined only by Xi we get that 

E(-^'-^^[exp{-^Ai{THmX2{TH))] = ii;(-i'-2)[exp(-^ H {X,{s) + X2{s))dsmX2{TH))] 



E-^exp{-^ f X2{s)ds)<P{X2{t)) 
^ Jo 



t = TH 



\2 l-TH 

^-i[exp(-— / Xi{s)ds) / ^(rH,X2,r)</.(r)dr] 



00 /'OO 



JO 



tt;(t, xi)V'(t, 2:25 r)(j){r)drdt, 



where 



and 



il){t,x,r) 



w{t, x) 



exp (-y ^ X2{s)d^ ;X2(t) G dr 



A^ rH 

exp( — — / Xi{s) ds) ■,TH £dt 
^ Jo 



From Lemma 4.2 putting u = —a/2, we obtain 



ip{t, X, r) = x"/^ 



) coth(iA) 



2sinh(a) 



xrA \ 



Vsinh(a) ) 



To evaluate w{t,x) we apply the formula (27) from Lemma 4.3 putting again v = —a/2. We obtain 



a2 



-L V 9 /I ~r 



£;-[exp(-7r^, - - Xi(s) d.)] = ^i;,l,/4r(»/2) ^"^'"'^'^-^/^^-"/^(^^^ 



(30) 



Now we invert the formula (30) with respect to 7. Using (|9| and the formula 25 p. 651 from [F] we obtain 
that 



w{t,x) = E"" 



X2 j-TH 

exp(— — / Xi{s) ds) ; T £ dt 
2 Jo 



^0/2^1+0/2 



2°/2r(a/2)sinhi+°/2(iA) 



exp 



xX cos 



h(U)\ 



2sinh(tA) J 



dt. 
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Combining all the above we get 



Jo 



w{t, xi)7p{t, X2, r)dt 



a a ^ 



2^+1 Til) Jo sinh2+t(a) 
After substituting A coth(tA) = 2s we get 



r""/'' / Axi + X2 + r) coth(tA) A ^ 
exp ^-A J I. 



sinh(tA) J 



(X2)4(X1)S 



-{xi+X2+r). 



22r(f )W4 A/2 

Now taking 2A instead A we have 



A2 



(X2)^(xi)t 

r(f)r°/4 



Coming back to initial variables (2/1,^2) = (4xiX2,X2 — xi) and (2/1,2/2) = (-Zi,-22) we obtain 



xi+X2 = Jyi + yi = \z\, xi = (|z| - Z2)/2 and X2 = (1^1 + Z2)/2 . 



Finally 



^(21,22) 



e 2 



(|z| +Z2)4(|z| -^2)2 

2^r(f)W4 A 

From ([5| we get that 



-{\z\+r)s( 



lim z'^I-i,(xz) = — ^ — . 

z^o "^^ ' r(l - u)x^ 



Consequently, for zi = 0, Z2 = n < we get 



e-^^^i;B{TD,) G dr 



(-2n)f2"/2 



2-r(f )r(i - |)2"/4W2 Jx 



sin(7ra/2) f —u 



"/2 p~X{r-u) 



r — u 



□ 



4.2 Hitting distribution of two half-lines in 

In this section we are interested in finding the distribution of the first hitting time and place of two 
half-lines for the process {r[~"^'^\ B{t)). 
Define the following two sets 

G = {(xi,X2) e [-1, 1] X [l,cx)) : -1 < xi < 1}, 
Ci = {(2/1,2/2) GMx [0,00) : (2/1 =0) A (|y2| > 1)}^ 
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and denote by tc^ the first exit time of the process {r[ °'^'^\B{t)) from the set Ci, i.e. tc^ = inf{t > 
: {rI~°'^'^\ B{t)) ^ Ci}. As previously we can consider the processes Y = (11,12) defined in (21 1 with 
Yi = (i?(""/2))2 Y2 = B. We thus have 

Tc, = inf{t > : Y{t) id}. 
Denote by tq the first exit time of the process X = (Xi,X2) with independent BESQ^~°'I'^'> components, 



as defined by (22 1, from the set G, i.e. 

TG = inf{t > : X{t) i G] = \ni{t > : |Xi(t)| = 1}. 
The function h was defined by hi{xi,X2) = (1 — xf){xl — 1). 




Figure 2: Simulated paths of X = {Xi,X2) and image Y = (li,l2) under mapping h. 



Using the same arguments as in Lemma 4.1 we get the equality between distributions of {tc\,Y{tc\)) 
and {A2{tg), h{X{TG))), where in this case the integral functional A2{t) is given by 



Mt) 



(Xlis) - Xf{s))ds 
(qiiXiis)) + q2{X2{s)))ds, 



where qi{x) = 1 — and q2{x) = — 1. Notice that qi{Xi{s)) ^ and 92(^2(5)) ^ almost surely. 

Lemma 4.5. The distribution of (rc^, y(rcj)) with respect to P^y^'y^^ is the same as the distribution of 

{A2{tg) , h(X (tg))) with respect to (xi,X2), where h{xi,X2) = (2/1,2/2)- 

Now we state the main theorem of the section. We give a representation of the density of {tg^ , y(rcj) 
in terms of spheroidal wave functions (see (36 1). 

Theorem 4.6. For {rI)~°^^\ B{0)) = (21,^2) £ Gi and r > X2 I we have 



e 2 ^c!i;5(rci) G dr 



1 



(r2 - l)"/2 2TTi J^/^_i^ 



m.^,x{xi)wxmi^{x2)(l)i ;,(r)di9, (31) 



where x, = \ ( V^f + (^2 + 1)^ + + (^2 - 1)^) , X2 = \ ( V^i + (^2 + 1)^ - ^ zl + {z2 - 1?) 
the function m^^x{-) is the solution of the following differential equation 



and 



(1 - x^)y"(x) - (2 - a)xy'(x) - (A^(l - x^) + 2'd)y{x) 



\x\ < 1, 



(32) 
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with boundary conditions m^^\{—l) = 0, m^^x{l) = 1. The functions 0^ ■^{■), (f>^ ^(•) are respectively 
increasing and decreasing independent positive solutions of the differential equation 



(r^ - l)y"{r) + (2 - a)ry'{r) - {X^^r^ - 1) + 2^)y{r) 
satisfying lima;^i+ = ^' ^^^x^oo 4>i^\{x) = and 



0, r > 1. 



(33) 



Proof. Let (f> G C^(l,oo) and A > 0. According to Lemma [45] we have 
We define the following functions 



(34) 



and 



^x{t,x,r) 



Wxit.x) 



exp(- 



A 



2 



q2{X2{s))ds)-X2{t)£dr 



, X ^ l,t ^ 0,r > 1, 



X2 rra 

exp(-— / qi{Xi{s))ds)]Xi{TG) = 1; tg G dt 
^ Jo 



\x\ ^ l,t > 0. 



Moreover, for every > 0, we define its Laplace transforms with respect to the variable t by 



M^,x{x,r) 



e-^^Wx{t,x)dt, 



e~^^^x{t,x,r) dt. 



Observe that 



e-'^^Wx{t,x) dt 





E^ 



A^ /-^G 

Xi(tg) = l;exp(-T9TG)exp(-— / gi{Xi{s))ds 

^ Jo 

Xi(TG) = l;exp(- r gl{X,{s))ds) 
Jo 



where ql{x) = ^(1 — x^) + t?. Using the Schrodinger equation (see |S] Theorem 9.10) we get that 

iGi'm^,x){x) - ql{x)m^^x{x) = 
{l-x^)ml^{x)-{2-a)xm'^x{x)-{>^^{^-x^) + 2'd)mi,,x{x) = 0. 

For x = —1 and x = 1 we have tg = and consequently ;^(— 1) = and m^^xi^) = 1- 
We have 



M^,xix,r) 



Observe that the function 



e'^^^xit,x,r)dt 



exp(- 



A 



2 rt 



q2{X2{s))ds)-X2{t)edr 



p{t,x,r) = E"" 



\2 ft 



exp(- 



q2{X2{s))dsy, X2{t) G 
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t > 0; X, r ^ 1, 



is the transition density of one-dimensional diffusion with the generator 

1,0 d?'^ 2 — a du 2 

From the general theory (see [S], Ch. II) we have that the T?-Green function of the diffusion is given by 

/•oo 

G^ix,r)= e-'''pit,x,r)dt = cxi^)-(l>l ^(x)^^^ ^{r), x<r, (35) 
Jo 

where ''^i a(^) positive solutions of the differential equation 

C2U — -du = 

such that (/)^ ^{x) is increasing and c/)^ -^{x) is decreasing and they satisfy the boundary conditions 
lima;^i+ 0,^ ^{x) = 0, lim.x^oo't'i a(^) ~ '^^^ Wronskian of the pair (0^ ^{x),(j)^ xi^)) given by 
= ^l,xi^)£(t^i,xi^) - £^l,xi^)(l^i,\i^) and the function cx{^) is given by 



cxi^) = wx{jl})fn{r) 
s'{r) 



wxm 



1 



m{r). 



(^2 _ x)a/2' 

where the function 

2 

wxi'^) = ? i 

does not depend on r. 

Recall that depends only on Xi and the processes Xi and are independent. Thus we easily 
obtain that the expression in the right hand side of (34) is equal to 

E-i[exp(-y j^J gi(Xi(s))d5)E^-nexp(-y ^*g2(X2(s))ds),/.(X2(t))],^=i,Xi(rG)=i] = 

/oo 
(t)[r)^x{t,X2,r)drdt 

/•oo 

0(r) / Wx{t,xi)^x{t,X2,r)dtdr. 



Observe that 

(•00 /•oo 



m^,x{x)\ ^ / e'^^^^^Wx{t,x)dt(^ / e-^W*Wo(t, x) = mjj(^)^o(^ 



and the integral defining ms({(^) o(a;) is finite for every ■d £ C such that 5R('!9) > —a/2 (see (38)). Conse- 
quently, using Fubini's theorem and the formula for the inverse Laplace transform we get 

/oo / 1 /■a/4+ioo \ 

Wxit,xi) f — / e'''M^^x{x2,r)d^]dt 

/•a/i+ioo / /•oo 

~ 27rz 



a/4— ioo V^O 



ci/4— ioo 

e^'Wx{t,xi)dt ) M^^x{x2,r)d^ 



ra/i+ioo 

m^^x{xi)M^^x{x2-,r) dd. 

a/A—ioo 



2m 
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Coming back to initial variables yi = (1 — xf)(x2 — 1), y2 = xiX2 and (yi,2/2) = {zf,Z2) we have 



(x2 - xif = yi + (2/2 - 1)^, {xi + X2)^ = yi + (y2 + 1)^ and X2 ^ xi. Finally 

^Zf + + 1)2 + ^Zf + (Z2 - 1)2) , 
/zf + (z2 + l)2-^z2 + (^^_l)2j 



and 



e 2 ^c!i;5(rcj G dr 



1 

27ri 



a/4+joo 



a/4— too 



m^i)^x{xi)M^^x{x2,r)dd 



1 



(r2 - l)"/2 27ri 7„/4_i^ 



a/4+ioo 



m_^,A(a:;i) • w\{-d)<j)\^^{x2)(t)'^^^xi.'^)d'&- 



□ 



Remark 1. The result in the case X2 > r ^ 1 can be obtained from the one given above by 
interchanging the role of X2 and r in the integral on the right-hand side of (31 ). This is an easy consequence 
of the symmetry of the Green function (35) with respect to the speed measure m{dr). However, the result 
given in (31 1 includes the most important case X2 = 1 which corresponds to zi = 0. 

e spheroi 

+ - z^Mz) = 



Remark 2. The equations (32) and (33) can be reduced to the spheroidal wave equation 

2 



(1 - z')y"iz) - 2zy{z) + [KAl) 



(1-Z2) 



(36) 



The radial spheroidal functions si]}i{-j,z) and sj^)i{'y,z) and angular spheroidal functions PS^^^{'y, z) and 
QSv^jiil^z) are solutions to the spheroidal wave equation in appropriate regions (see [21J). When 7^ = 
the equation (36) reduces to the Legendre equations (14). When = 1/2 (i.e. a = 1) the equation (36) 
can be reduced to the Mathieu equation and the spheroidal wave functions are reduced to the Mathieu 
functions (see [TH], [T]). 

To illustrate our method we compute what appears to be the Poisson kernel of the interval [—1,1] for 
the standard isotropic a-stable process (see Proposition 3.1). 

Corollary 4.7. For A = 0, zi = 0, |z2| < 1 and r > 1 we have 



E(0'^2) [B{tc,) e dr] 



sin(7ra/2) f 1 — z. 



TT 



2\ a/2 



1 



Z2 



Proof. For A = the equation (32) becomes 

(1 - x^)y"{x) - (2 - a)xy'{x) - 2i?y(x) 
Making substitution y(x) = u{t) with t = {x + l)/2 we get 



0. 



t(l - t)u"{t) - [1 



a 



(2 - a)t]u{t) - 2^u{t) = 0, 



(37) 



i(l - a) + i((l - a)2 - 8i9)V2, 
dng into account the general solution of and the boundary 



1 — ^ and a 



which is the hypergeometric equation (12) with c 

h=\{l-a)- - a)2 - 8i?)V2. Ta: _ 
conditions m^(— 1) = and m^(l) = 1 together with (11) we get the following formula for q 

'l + x\"'^T{a/2 + A{'d))T{l + a/2- A{^)) ^ , , , , ...... ,^l + x. 



m^{x) 



r(a/2)r(l + a/2) 



2Fi{A{^),l-A{^y,l + a/2; 



(38) 
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where 



Substituting y{r) = (r^ — l)°'/^w{r) in (33) with A = we arrive at the Legendre differential equation 
(14) with /i = a/2 and v = B{-d) = — | + ^((l — a)^ + S-!?)^/^. Consequently, the monotone solutions can 



be chosen as 



4oW = {r^-^r"Pl[l){r), 



From ( 15 ) we get that 



'^l.o(l) 



2^/2 



r(l - a/2) 



and 



2^/2 
r(l - a/2)' 

To evaluate the integral appearing in the right-hand side of ([31]) we integrate the function 

d m_^{xi)wo{^)(t)l Q{l)(l)\^Q{r) 
over the following contour of integration 



(39) 




The function of complex variable ■(? — > m_^(a;) is meromorphic in the half-space K(??) > with poles at 
points ■!?„, such that 

1 + a/2 - A{-^n) = -n 
what gives ??„ = |(a -I- n)(l -|- n). We have 



Res^„(r(l + a/2-^(-i?))) 



lim (i? - 7?„)r(l + a/2 - A{-^)) 

nm 

r(-a/2 ^(-^?)) sin7r(l a/2 - A(-^?)) 
.n l + a + 2n 
^ ' 2r(n + l) ■ 
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Consequently 

Res ^„ (m_^(x) wo(??)(/'l^ o(l)'^i, = 

-l)"r(l + g + n) + f + n, -n - f ; 1 + f ; . , 
2r(n + l)r(a/2)r(l + a/2)(l + a + 2n)-i ^o(^n).^,„,o(l)'^.„.o(0. 

Moreover we have ([E] p. 105 2.9.(18)) 

l + x\°^^ ^, , , , 1 + x, 



2Fi(l + a/2 + n, -n - a/2; 1 + a/2; 

1 — x\ ^, ,1 + X, 



2Fi(-n, n + 1 + a; 1 + a/2; 
r(n + l)r(a + l) 



r(n + l + a) V 4 

(-l)"C(^)(x), 



r(n + l)r(a + l) 



r(n + l + a) 

where p = (1 + a)/2 and C^^ is the Gegenbauer polynomial. We also have B{'dn) = n + a/2 and the 
Wronskian of the pair {4>\^ q, 0:^^ q) is given by 



(^2 _ -|^^a/2-lgja7r/22a" ^ 2 



r(^^ + i)r(^^±f±i) 
r(i + §)r(^^ 



Thus 



(2 .^a/2-l W2 r(n + a+l) 

^ ^ r(n + i) ■ 



2 

^O(^n)</'l„,o(l)4„,o(0 = ,2 n»/2-lM/J^T , ,, 4„.o(l)'/'l,o(0 

= ,-W2 r(n + l) 2"/^(r^-l)"/^ /2 . . 

r(n + a + l) r(l-a/2) '°^n+a/2V )■ 



Combining all above we conclude that the residuum of the function (39) at point 'dn is equal to 

-W2 r(a + l)(l + a + 2n) / l-xf y^' r(n + 1) , _ ^^a/4^«/2 . x 

r(a/2)r(l + a/2)r(l-a/2) V 2 ^ ^^^J ^ ^ ^ ir ij ^^^^/^yr) 

- e-W2^«?!^(l _ .?)"/^C(^)(x2)^(^^(.^ - ir^^Q-J^r). 

Letting R ^ oo and using the asymptotic expansions of the considered functions we obtain that the 
integral of the function (39) over the line (a/4 — icx),a/4 + ice) is equal to the sum of the residues at 
points -dn- Details are left to the reader. Consequently, using (31 1 we get that E'^^'^^) [b{tci) S dr] is 
equal to 

sin(;ra/2) (1 - zj)"/^ ^ (1 + « + 2n)2tr(i±^) r(n + l) , _ /, . .^(,). . 

7, (^2_l)a/2Z^^ ^ r(n + a + l)^ ' ^n+a/2^^)^n i^2j 

The relation (P] 7.312 (1)) 

"1 (l_^2)a/2 ^2W2 _^ ,,a/4^a/2 , , 
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and the orthogonal relations for the Gegenbauer polynomials 



give 



where 



nT{n + 1 + a) 



-1 (n + ^)r^(^)r(n + 1) 



^ oo 

— = ^a„(r)C(^)(x), 



r — X 

n=0 



^n[r} - e ^^^^ (r Ij Q^^^/,[r} 2 ^r(n + 1 + a) 

_ ^ (2n + a + l)2"/M^) Tjn + l) 



Xn + l + a) 
Using the above-mentioned formula we get 

sin(7ra/2) (1 - z2)"/2 i 



E(0'^2) [^(^^j G dr] 



TT (r2 - l)"/2 r - Z2' 

□ 



5 Hitting distributions in 

The aim of this section is to generalize the results of Section 4 to higher dimensions. Assume that < a < 
2. Let Y{t) = {Yi{t),B"'{t)) be the Brownian-Bessel diffusion in M"+i defined in Section|3| That is Yi{t) is 
a BES''~"/'^^ process independent of n-dimensional Brownian motion B'^{t) = {B2{t) , B^{t) , . . . , Bn+i{t)). 

We consider the (n-l-l)-dimensional process Y{t) exiting from various open sets which are complements 
of lower dimensional subsets of R"+^. We define 



Dn 


= {yGl 


: (2/1 


= 0) A {y2 > 0)}^ 


Cn 


= {yGl 


: (yi 


= 0)A(|y2| >l)r 


Hn 


= {yGl 


: (2/1 


= 0)A(y2 = 0)r. 



Observe that the above sets are complements of an n-dimensional halfspace, an n-dimensional strip and 
an (n — l)-dimensional linear subspace, respectively. 

Throughout the whole section we use the following notation. For a point x = (xi, 3:2, . . . , Xn+i) € M"+-'^ 
we denote hy x = (a;2, . . . ,Xn+i) its projection onto M", and by x = (xs, . . . its projection onto 

We also denote points from M" by x, y, etc. Likewise points from M"~^ are denoted by x, y, etc. 
For any points y = (yi, 2/2, • • • , Vn+i) G M"+\ x = (x2, . . . , x„+i) G and x = (xs, . . . , x„+i) G M""^ we 
denote by 

/ n+l \ / n-l-1 \ 

|y - x| = I + ^ |yj - I , |y - x| = [yl + yl + ^\yi- xi\^ 1 

the Euclidean distance between the point y G M"+i and (0,x) G M"+i and y G M"+^ and (0,0, x) G M"+i, 
respectively. 

We define the first exit time of Y from Dn by 

rz)„ = inf{t > : Y(t) ^ 

We begin with providing the formula for the joint density of (r/j^, Y(t£)^)) when Y starts from the 
point y = (yi, . . . , y„+i) such that yi = 0. 
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Theorem 5.1. For y = {0,y2, ■ ■ ■ , Vn+i) such that 1/2 < we have 

mro.. . nr.j . M ^ « (?)°'^ ii^»p {-'^f) . («) 

where t > Q and a = (cj2, . . . , (Jn+i) G , o"2 > 0. 

Proof. We begin with the case n = 1. Then applying the following formula for the Laplace inverse 
transform (see [5] formula 2. page 650) 

a f a?'\ 



to the formula (29) with A replaced by -y/27 we obtain 

fy'{t,a2) = E^^'y^^ [td, e dt,YiTD,) G da2] 
= E^^'y^^ [td, e dt, W2{td,) G da2] 

sin(W2) ("-^^r^'^-i / e-v^(-^-^^) \ 

sin(7ra/2) /-y2V^^ ^ (0-2-^2)^ 



2V2vr3/2 ^2 ; t3/2 \^ 2t 

This ends the proof for n = 1. When n ^ 2 we observe that the density of the joint distribution of 
{td^, B2{td„)) with respect to P^^^v^) is equal to /^^(t, (T2). Let gt{x) denote the density function of the 
distribution of (n — l)-dimensional Brownian motion. Obviously, we have 

9t{x) = exp — — 

{2TTt)— \ 2t 

Using the fact that the first exit time td„ depends only on Yi and B2 we obtain 

Ey[TD„ G dt, Y{tdJ G da] = Ey[TD„ G dt, B2{tdJ G da2, BsirnJ G das, ■ ■ ■ S„+i(rD,J G dan+i] 

= E^'^'y^^im^ G dt, B2{tdJ G da2] gt{y - a) 
= fy\t,a2)gt{y-a) 

sin(7ra/2) ( -y2\"^'^ 1 ( l^'-yP^ 

□ 



2n/2^1+n/2 \ a2 I t^+^l"^ V 2t 



Now we present the multi-dimensional generalization of result given in Theorem 4.4 
Theorem 5.2. For y = M""^-*^, such that (1/1,1/2, • • • , J/n+i) G -Dn we have 

[e ,rirBjGflaj (27r)"/22°/2r(a/2) 

4sin(W2) a,n+a/2 f / "^2 \ ^ ^=±^ (^1 ^ " ^n/2( A|cr " 5| ) 
yi'^ / 1 I ^ r: ::^77r7^ 



2nW2^n+lr(«/2) «i 7(_oo,0)xM-l \ ^2 J |y - 1^ " ^l"/' 

where a = {a2, ■ ■ ■ , o"n+i) G M", (T2 > 0. For yi = we get 

2si n(7ra/2)A"/2 / -^2 V^' ^n/2(A|y - a |) 
2!^=^" U2y 



YimJ G da] = ----:^r(^^ ( ^ ) "^^ 'l..; ■ (41) 
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Proof. For y = (yi, 1/2, • • • , Vn+i) G Dn such that yi = we use Theorem 5.1 and formula ^ to show that 



sm(7ra/2) ( -y^y'^ .2 



n n-\-2 

227r 2 



<T2 / Jo 

2 sm(7ra/2) ,^ f -y^ \ i^n/2(A|a - y| ) 

For the general point y = (yi, y2i • • • 1 yn+i) = (2/1, y) G M"^-'^ we define the first hitting time of by 
the process Yi, 

T = mi{t > : Yi{t) = 0}. 

We denote by h{yi,t) the density of T. This density was found by Getoor and Sharpe (see [12]) and is 
given by 



h{yi,t) 



Vi 



exp 



yl 



Using the independence of T and and ^ we easily obtain 



t > 0. 



e 2 



■T. nn 



e 2 



* ^ exp(-^— 



(27rt)"/2 



2t 



e 2 



(2^)"/22"/2r(a/2) io t^Hn+a)/2 



exp( — )dt 



2yfA- 



-R'n+Q (A|y — z\) 



(2^)"/22"/2r(a/2) |y-z|^ 
Using the strong Markov Property we obtain 



Ey 



e 2 

Ey 
Ey 



TDr, 



;Y(rDj e da 



+ Ey 



T = TD„;e-^^^«;Y{TDj €da 



e-^-On;Y(TDj G da 



+ Ey 



B2{T) ^ 0;e-^^;Y(r) G da 



Consequently we get the following expression for £'^[e 2 "^^n; Y(r£i,J G da] for every y G 

/ Px{y,2)Hx{z,a)dz + Px{y,a). 

J(-oo,0]xR"-i 

This ends the proof. □ 



The following corollary is an immediate consequence of Theorem 5.2 (formula (41)) and Proposition 



Corollary 5.3. Let Dn C M" be the halfspace {x G i?"; xi < 0} C M". T/ien t/ie m-Poisson kernel of 
Dn for the relativistic a-stable process with parameter m > is given by 



2sin(7ra/2)m2"« "/2 Kn/2{m^ \y - o- 



n n-f-2 

2 2 7r 2 



0-1 



\y-o- 



n/2 



20 



where y = (yi, . . . , y„) € -D„ and a = (ui, . . . , cr„) G Z)^. For m = we obtain the Poisson kernel of Dn 
for the standard isotropic a-stable process given by the formula 



sin(7ra/2)r(f) 



n + 2 

vr 2 



-yi 

0-1 



a/2 



1 



\y - 



The conclusion of the above corollary (for m > 0) is one the main results obtained in [B], by using a 
different approach, which was almost entirely analytical. The present presentation provides another proof 
which is of probabilistic nature. For m = 0, i.e. for the standard isotropic a-stable process, the formula 
for its Poisson kernel of Dn is well known and can be obtained from a formula of the Poisson kernel for a 
ball (see P]), where an indispensable tool was Kelvin's transform. 

We define the first exit time of Y from C„ by 

rc„ = inf{t > : Y{t) i C,}. 
In order to describe the distribution of (tc„, Y(tc„)) we define 



e 2 ^Ci - ^2(rci) G dr 



\r\ > 1, 



see 



where {zi,Z2) G Ci. The integral representation of H{zi, Z2, \,r) is the main result of Section 4.2 
the formula ([31])). We define 

/i(A,y,a) =^^[e-^"C";B"(rcJ Gda], a = (da, . . . , G M", l^al > 1. 

In the next theorem we provide a formula for the (n — l)-dimensional Fourier transform of /i(A, y,o") 
which entirely describes the distribution of (rc^, Y(rc^)). 

Theorem 5.4. Let n > 2. For y = {0,y2, ■ ■ ■ ,yn+i) such that \y2\ < 1 and z G M""-*^ we have 

[ hiX, y, a)e'^^''^da = i/(0, ya, V|y-^P + A2, ^2). 

Here a = (o"2, cj). 



(42) 



Proof Let fy^{t,a2) = E^'^'^^'i [tc^ G dt,B2{Tc^) G da2]. When n ^ 2 we observe that the density of the 
joint distribution of {tc„ , B2{tc„)) with respect to p(0'2^2) equal to /^^(t,cr2)- Let gtix) denote the 
density function of the distribution of (n — l)-dimensional Brownian motion. Using the fact that the first 
exit time tc„ depends only on Yi and B2 we obtain 

Ey[Tc^ G dt, Y{tcJ G da] = Ey[Tc^ G dt, B2{tcJ G da2, Bs{tcJ G das, ■ ■ ■ Bn+i{TcJ G da^+i] 

= [tc„ G dt, B2{tcJ G da2] gt{y - a) 

= fy\t,a2)gt{y-a). 

This implies that 

a2 a2 

h{X,y,a) = Ey[e-^-c^;Y{rc^)^da]= / e-^'fy^{t,a2)gt{y-a)dt. 

Jo 

Hence, the formula for the (n — l)-dimensional Fourier transform of h{X,y,a) can easily be found 

/ /i(A,y,a)e^(^'^")(ia= / e-^V^^(t, ^2) / gt{y - a)e'^~^''^dadt = / /^^(t, (J2) e-*^'^^2 dt. 

Jk"-i jo jir"-i jo 

Next, observe that the last integral is equal to H{0, y2, y/\y — + A^, 0-2). □ 
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Again, the following corollary is an immediate consequence of Theorem 5.4 and Proposition 3.1 



Corollary 5.5. Assume that n > 2. Let Cn C M" he the strip {x e M"; |xi| < 1} C M". Then the 
(n — \)- dimensional Fourier transform of m-Poisson kernel of Cn for the relativistic a-stable process with 
parameter m > is given by 

where y = {yi, ...,yn)^Cn and a = {ai,a) £ C^. 

Note that the above formula for m = provides the (n — l)-dimensional Fourier transform of the 
Poisson kernel of the strip Cn for the standard isotropic o-stable process. 

In the last part of this section we assume that 1 < a < 2. We consider a similar problem as above 
but now we compute hitting distributions when the (n + l)-dimensional process hits (n — l)-dimensional 
half-space, with n > 1. Let 



{y G 



p71+l 



: (yi = 0) A (y2 = 0) A (ys > 0)}^ 



and let th^ be the first exit time of Y from Hn- 

We will reduce the problem to an n-dimensional situation. Observe that Z = \/Y1 + B\ is a Bessel 
process of index (1 — a)/2 (see [2DJ, Ch. XI) so W = {Z, B-^, . . . , Bn+i) is our n-dimensional Bessel- 
Brownian diffusion. Obviously all of its components are independent. Let 



Hri 



{y G M" : {yi = 0) A (ya > 0)}^ 



We define Tjj^ being the first exit time of W from Hn- Observe that Y exits Hn if W exits Hn- Moreover 

(r^^,(0,W(r^J)) = (TH„,Y(r^,J). 
Hence our problem is within the framework of the problem studied in the first part of this section (for 



the process W) but with n replaced by n — 1 and a by a — 1. Applying Theorem 5.2 (Theorem 4.4 if 
n = 2) we obtain 



Theorem 5.6. For y 



nn+l 



, such that yf + y2>0 and 1 < a < 2 we have 



Ey[e-^^"«; Y{thJ G da] 



2(^2 ^ y2^^^^ K^{X\y - a\) 

n—1 a — 1 , „ 1 , , , n — 2 + q: 

{2TT)^2^r{^) \y-a\^^ 



+ 



+ 



4sm{7T^)iyi+yi) — X 



a~l 2n--~3-{-a 



-Z3 



Kr 



,{\\y-~z\)Kr^{\\a--z\) 



\y-z\ 2 



-I 

a — z\ 2 



2 2 7rnr(^) J(-oo,0)xM"-2 V CIS 

where a = (0-3, . . . , ct„+i) G M"~^, 0-3 > 0. For yi = y2 = we get 

2sin(^(a-l)/2)A("-i)/2 ^-yaA ^(n-i)/2(A|y - ct|) 



Ey[e-^^""; Y{thJ G da] 



n — 1 n-fl 

2 2 vr 2 



\y _ o-|("-l)/2 



For n = 2 the first formula can he simplified to 



e 2 '"2 



; B3{th2) G dr 



/II \ /II \ 

{\z\+Z3) 4 (l^l - Z3) 2 



where \z\ = \J z\^ z\^ z\ and z\^ z\> 0. 



-dz. 



g-(|.|+r)s(^2 _ ^2)^^^ (^^27^1^1 +^3\/s2 - ds, (43) 
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Again, the following corollary is an immediate consequence of Theorem 5.6 and Proposition 3.1 



Corollary 5.7. Let H2 C be the complement of the half-line {x £ M?; xi = 0, 3:2 > 0} C M^. Then 
the m-Poisson kernel of H2 for the relativistic a-stable process with parameter m > and 1 < a < 2 is 
given by 



{\y\ +y2) 4 {\y\ -y2) 2 
2^^r(^^)r("-i)/4 Jmi 



<\y\+r)s(^2 



where r > and y = (yi, 2/2) S H2 with yi / 0. 
If yi = and y2 < we have 



Pgiy,r) 



sin(7r(a - l)/2) 



-y2 \ ^"""^^/^ g-ml/Q(r-M) 



r — u 



For m = we obtain the Poisson kernel of H2 for the standard isotropic a-stable process given by the 
formula 



ce — 1,, , ,CK — 1 



2-4-^r(^^fi)r(°-i)/4 Jo ^ ^ ' 

where r > and y = (yi, 1/2) G -^2 wit/i yi / 0. 
If yi = and y2 < we have 

2 TT \ ^ / r — u 

To our best knowledge the above formulas have not been known before even in the stable case. As far 
as we know the only explicit formula related to the standard isotropic a-stable process (1 < a < 2) killed 
on exiting H2 was the formula for the Martin kernel of H2 with the pole at infinity obtained in [J]. 



Next, we present a multidimensional version of Corollary 5.7 For the purpose of clarity we give the 



description of the Poisson kernels if the starting point of the process belongs to the subspace spanned by 



the underlying half-space. The general case can be easily recovered from Theorem 5.6 



Corollary 5.8. Let Hn C M" be the complement of the (n — 1)- dimensional half-space {x G i?"; xi = 
0,X2 > 0} C M". Then the m-Poisson kernel of Hn for the relativistic a-stable process with parameter 
m > and 1 < a <2 is given by 



p'S (y^ ^) 



2sin(7r(a-l)/2)m5s f -y2Y" ^^^^ i^(n-i)/2("T'" |y 



2^7r^ |y_^|(n-l)/2 



where y = {yi, . . .,yn);yi = 0,y2 < and a = {(T2, ■ ■ ■,crn); (72 > 0. 

For m = we obtain the Poisson kernel of Hn for the standard isotropic a-stable process given by the 
formula 

sm(j(a-l)/2)r(fe^) f-y2\"-'>'^- 1 



\y — cj|" ^ 



It is very striking that the above formulas are identical with the formulas of Poisson kernels of {n — 1)- 
dimensional half-spaces for the stable process or relativistic stable with index a — 1 living in M"~^. But we 
have to keep in mind that this is true only when the n-dimensional process (either a-stable or relativistic 
a-stable) starts from the subspace spanned by the complement to Hn- If the process starts from other 



points the formulas become much more complicated (see e.g. (44 1 for the general two-dimensional stable 
case). 
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6 Appendix 



6.1 Proof of Proposition |3.1 



Recall that Y(t) = (Yi(t), is the Bessel-Brownian diffusion in 11"-+^ defined in Sectionjsj We begin 

with computing the A-resolvent kernel of the semigroup generated by Y{t). Recall that the transition 
density (with respect to the speed measure m{dy) = 2y^~°'dy, y > ) of Yi{t) is given by the following 
formulas: 



and 



(-a/2) 



(0, y) = 2"/2-it°/2-ir(l - e 



For a = 1, Yi = ^/2) ig ^j^g reflected Brownian motion with the density 



Qt 



(-1/2) 



^-\x~y\^/2t _^ ^-\x+y\^/2t 



(45) 
(46) 

(47) 



We use the following notation: y = {yi,y), where y = (?/2, • • • ,yn+i) G M". Hence we find that the 
transition density of Y{t) with respect to the measure m{dyi) x dy,yi G M+,y G is equal to 

q{t;x,y)=ql-''/^\xi,yi)gt{x-y), x, y G M'^+S 

where (x — y) is the transition density of -B". Observe that q{t; x, y) is symmetric. Let C/a , A > denote 
the A2/2-resolvent kernel for the process Y(t). That is 

/>oo 

f/A(x,y) = / e'^^*/^q{t;x,y)dt. 
Jo 

Lemma 6.1. Suppose that x = (0,x) or y = (0,y) . Then 

(A/2)^ K^{X\x-y\) 



Ux{x,y) 



W2r(l - a/2) 
For a = 1 we obtain, for all values of xi,yi: 



(48) 



\x — y\ 2 



Uxix,y) 



1 /^\ 2 
vr V 27r 



Kn^{X\x-y\) Kn^{X\x-y*\) 



\x — y\ 2 



+ 



n-1 

|a; — y*! 2 



(49) 



where for y = (yi, y) we have y* = (— yi, y). 



Proof. By symmetry of C/a we may assume xi = 0. Using ( |19[ ) we obtain 

(2t)°/2-i 



Ux{x,y) 



e 2t 



r(l - a/2) (2vrt) 
(2t)°/2-i 



n/2 
2 



2t e 2 (it 



e 2t A^i 

e 2 est 



r(l - a/2) (2vrt)"/2 

(A/2)'^ i^:i^(A|x-y|) 



W2r(l - a/2) 



Ti — a 

\x — y\ 



For the case a = 1 we use the formula (47) to obtain (49 1. 



□ 
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Now, let D be an open set of M""*"^ and td = inf{t > 0; Xt ^ D}. A point y € D'^ is called regular for 
D" if 

PV{TD = 0) = 1. 

Note that by Blumenthal's — 1 law the equivalent condition for regularity is -P^(tz) = 0) > 0. It is clear 
that all points in the interior of D'^ must be regular. Hence regularity must be justified only for points 
from dD. From the general theory it follows that if y € is regular for D and x S D then 

Ux{x, y) = E^[td < oo; e~'^^^U^{Y{TD),y)]. (50) 

Next, we consider sets whose complements are of lower dimensions. Let F be a closed subset of M" 
and take D = M"+i \ {0} x F. Observe that 30 = 0" = {0} x F. We want to impose a condition on the 
set F which guarantees that a point of dD is regular. We say that F C has the interior cone property 
at y G F if there is an open cone F C M'^ with vertex at y and r > such that Fr = T n B{y^ r) C F. 
Note that in the case n = 1 the above property is equivalent to the condition that for a point y G F there 
is an open interval ly with one of the endpoints equal to y such that ly C F. In particular if we take 
F = /„, where I„ are disjoint closed intervals, then every point of F has the above property. 

Lemma 6.2. Let y = (0, y) G {0} x F and F has the interior cone property at y £ F. Then the point y 
is regular for D'^. 

Proof. We may and do assume that y = 0. Let % = inf{t > 0, Yi{t) = 1/k}, = mf{t > Yi{t) = 0}, 
where A; = 0, 1, We claim that 

P°( lim Tfc = 0) = 1. (51) 

Let e > 0, then from the scaling property for the Bessel process we have 

pO(Tfc > e) = P°(ri > fc^e). 

By the fact that Yi hits any point of [0,oo) with probablity one we infer that P°(Ti < cc) = 1, which 
implies that the decreasing sequence is convergent in probability, hence P° a.s. Observe that on 
the set {Y{Tk) G {0} x F} = {B''{Tk) G F} we have td < Tk. Hence from ^ we infer that 

F°(rD = 0) > P°(limsup{S"(rfc) G F}) > limsupP°(P"(rfe) G F). 

k k 

Next, 

p\B'''{Tk) eFnr)> p\B^{Tk) g r^.r^ < i). 

Noting that P°{B"{t) G T^) > F°(P"(1) G T^) > 0, < t < 1, and using independence of the process 
B"^ and we obtain 

lim sup P^iB"^ in) G F) > P^{B^ G Tr)P^{Tk < 1) ^ F°(S^ G T^) > 0. 

k 

This implies that F°(tz) = 0) > and then P'^{td = 0) = 1. □ 

Let X"^ be the a-stable relativistic stable process living in M". The kernel of the A-resolvent of the 
semigroup generated by X'^{t) will be denoted by U^{x,y). We have 

;>oo 

U^{x,y)= / e-^'pT{x-y)dt, 
Jo 

where Pt^{x — y) is the transition density (with respect to the Lebesgue measure) of the process X™ . 
The function C/™ has a particularly simple expression when A = m (see |6j): 

rr^(. ,A - 2^-('^+")/^ K^^.^y,{my-\x - y\) 
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For open D C M" we define = inf{t > 0; X"^{t) ^ D} . From the general theory it follows that if 
y £ D'^ is regular for D'^ (for the process X™) and x G D then 



Recall that the m-harmonic measure of the set D is defined as 



for a Borel set A. If all points of D'^ are regular then the equation (53) can be rewritten as 



/ U:::{z-y)P^{x,dz) = U:;:{x-y), x^D, y^b^ 

J D'^ 



(53) 



(54) 



(55) 



which together with the following uniqueness lemma (see |S]), is crucial for obtaining the explicit form of 
P^{x,dz) in the particular cases studied in the paper. 

Lemma 6.3 (Uniqueness). Suppose that /i is a finite signed measure concentrated on a closed set F C 
with the (finite energy) property : 




F JF 



K(n-a)/2{m^'°'\z - y\) 
\z — y|("-")/2 



\li\{dz) \^\{dy) < oo. 



If for every z £ F we have 



then /i = 0. 



\z-y 



(56) 



(57) 



Comparing (52) with (48) we see that for points x = {0,x),y = (0,y) G M"+^ we have the following 
equality: 

U;;^{x,y)=CaU^^/4x,y), (58) 

r(i— -)2i~" 
where Ca = . 

Lemma 6.4. Let D d he open. Assume that F = \ D has the interior cone property at every 
point. Let D = M"+i \ ({0} x F). Let x = (0,x) £ {0} x D . Assume that P''{td < oo) = 1. Then the 

measure P{x,da) = E^{e 2 '^^ ] B'^{td) S da) is the m-harmonic measure of D for the n-dimensional 
a-stable relativistic process with parameter m > 0. 



Proof. By Lemma 6.2 all points of D'^ = {0} x F are regular for D'^. Then by (50) the measure P{x,da) 



concentrated on F has to satisfy 

Umi/cix,y) = j U^i/c{a,y)P{x,da), y £ {0} x F. 



By (58) this is equvalent to 



U:^{x,y) = I U::{a,y)P{x,da), yeF 
Integrating both sides with respect to P{x, dy) we obtain 



(59) 



U:^{i,y)P{x,dy) 



Um{cF, y)P{x, da)P{x, dy). 



Since x £ D, then supy^pUmi^yy) < °° the integral on the left-hand side is finite. This implies 
that the sub-probability measure P{x, da) has the finite energy (defined in Lemma 6.3 ). Observe that the 
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set F C M"" has the interior cone property, therefore all points of F are regular (for the process X™). The 
proof of regularity is identical or even simpler than the proof of Lemma |6.2[ Instead of the sequence of 
random times Tk, k > 1, one can take a detrministic sequence 1/k and use the fact that the process X"^ 



is rotation invariant. This implies that (55) holds and by the same arguments as above, for the measure 
P{x,da), we infer that the measure P^{x,da) also has the finite energy. Finally we may apply Lemma 



6.3 



together with (|55|) and (|59|) to claim that P{x,da) = P^{x,da). 



Observe that the above Lemma provides the proof of Proposition |3.1| for m > 0. It remains to settle 
the case m = and this is done below. 

Corollary 6.5. Let X{t) be the standard isotropic a-stable process in M". With the assumptions of 
Lemma \6m we have 

P'{X{r^) G A) = P-{B^-{td) e A), 

for every Borel set A G M". 

Proof. It is known that we can represent the stable process as X{t) = X™-{t) + Z"^{t), where Z"^{t) 
is an independent of X"^{t) compound Poisson process with interarrivals of jumps having exponential 
distribution with parameter m (see [22]). This implies that there exists T = T{m) with exponential 
distribution with parameter m, independent of X"*, such that 



Let 



T 



X{t) = X"\t), t<T. 



r]{m) =mf{t > 0,X"'{t) ^ D} 
mi{t > 0,X{t) i D}. 



(60) 



Observe that (60) implies that for every Borel set A G M" we have 



{X™(r?) €A,r]<T} = {XrGA,T< T}. 
Then for every Borel set A we obtain 

\E^[e-'^v^X"'{i]) eA]- P^{X{t) G ^)| < |S^[e-"^ X™(r/) eA,rj<T]- P^{X{t) e A,t < T)\ 

+ P^{7] > T{m)) + P*(r > r(m)) 

= |S*[e-"^X(?7) e A,T] < T] - P^iX'^ir]) e A,r, < T)\ 
+ 2P^{r] > T) 
< 3S*[l-e-n- 

In the last step we u sed the independence of rj and T. By the relationship between the processes Y and 

X"^ (see Lemma 
this implies that 



6.4 



we have E^e = E^e 2 '^o ^ which yields limm^o-E'^[l — e "^^] = 0. Finally, 



P^{X{t) G A) 



lim E^[e-'^'^,X'^{rJ) G A] 

m— »0 

"^;i?"(TD) G A] 



= lim E^[e "'2 

m— »0 

= P"(i?"(rD) G ^) 



□ 
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6.2 Generalized Bessel processes 

In this section we consider the process X = (Xi, X2) given by the following stochastic differential equations 

dXi = - Xl\ dBi - ^-^^1 dt (61) 

2- a 



dX2 = ^J\Xi-l\dB2 + ^-\X2\dt (62) 

such that Xi{0) = xi, ^2(0) = X2, where ^ 1 and X2^ I. Here Bi, B2 are two independent standard 
Brownian motions in M. We claim that that for 1-^^1(0)1 ^ 1 we have |-'^i(t)| ^ 1 for all t > 0. Analogously, 
X2{t) ^ 1, whenever ^2(0) ^ 1. 

The first process is a "quadratic" version of Legendre process; the second one is "'quadratic" hyperbolic 
Bessel process (see [2U], Ch. VIII, p. 357). Indeed, changing variables Xi = sinli and X2 = coshy2 we 
obtain 

1 - a 

dYi = dBi ^tanyidt (63) 

1 - a 

dY2 = dB2 + ^—cothY2 dt (64) 



with the generators 



1 d^ 1-a d 
1 1 — a , d 



These processes were investigated by various authors, among them by [15] but for different values of 
coefficients appearing in the drift term. Therefore, we present here, for convenience of the reader, a more 
detailed information about behaviour of these diffusions. 

Observe that the processes Xi^X2 exist as solutions admitting explosions (see Theorem 2.3 in ^17j, p. 



159). Moreover, the corresponding functions ai and hi in the equations ([61| and (62 1 

a,{Xi{t)) dBiit) + bi{X,{t)) dt, i = 1, 2 

satisfy in both cases 

|c7i(2;)|2 + |6i(x)|2^2(l + |x|2), 

which shows that the explosion time C = 00. 

Now, we prove that for |Xi(0)| ^ 1 we have |Xi(t)| ^ 1 for all t > and the boundary points 1 
and (—1) are instantaneously reflecting. Analogously, X2{t) ^ 1, whenever ^2(0) ^ 1, and the boundary 
point 1 is also instantaneously reflecting. We only sketch the first part of the statement. To do this, 
observe that local times Lj{Xi), L^^(Xi) are equal 0: 



t-l> 



l(i<Xi{s)) ds = l(i<xi(s)) 1^ _ jj^(g)2|^ = \l-a'^\ (^1) '^'^ • 
Thus, Ll{Xi) = 0. Analogously, L^^{Xi) = 0. By Tanaka Formula (see [20], Ch. VI) we have 

(Xi(t) - 1)+ = (Xi(0) - 1)+ + l(i<xi(.)) dXr{s) + \l\{X^) . 
Applying L\{Xi) = we obtain whenever |Xi(0)| ^ 1 

(Xi(t)-l)+= f l(^^x,(s))dXi{s). 
Jo 
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Taking expectations we get 



E{xi{t) - ly 



a 



E[l 



(l<Xi(s)) 



Xi{s)]ds. 



Now, the left-hand side is non-negative while the right-hand one is non-positive so both are 0. Conse- 
quently, for every t > we have Xi{t) ^ 1 almost everywhere. The same arguments apply to show that 
Xi{t) ^ —1 almost everywhere. Thus, the absolute values in equations (61) and (62 1 can be discarded. 
Another application of local times shows that 



LliX^)-Ll-iX,) = -2- 



a 



L(Xi(s)=l) 



ds 



Since obviously = we obtain that that the time spent by Xi in 1 has zero Lebesgue measure, 

which shows that the point 1 is (instantaneously) reflecting. The same conclusion holds true for the point 
(— 1) and for the point 1 for the process X2. 

Now, with the information that |^i(t)| ^ 1 for all t > we can prove the uniqueness of the process 



Xi{t). To show this observe that IvT" 



— \/ i — ^ \/\x^ — y^l = y^\x + y\ y^x — y\, which means 
that the function ai is locally Holder continuous of exponent 1/2. Moreover, the function bi is obviously 
Lipschtiz continuous. This assures the uniqueness of the process Xi{t) (see remarks after Theorem 6.1, 
p. 201 in [T7j). The same observation is true for the process X2{t). 

Obviously the processes Xi, X2 are independent. Moreover, the generators of the processes are given 

by 



n 



2- a d 
~ir''di' 



2 

_x2-l(i2 2-a d 



(67) 
(68) 



respectively. We assume that the domain of (67) consists of the functions u G C'^[—\, 1] with the property 
that n'(— 1) = u'{l) = and, in the case of (68), we assume u G C^[l,cx)) along with the property 
u'{l) = 0. Derivatives here are meant as (appropriate) one-sided ones. 

To classify the boundary points we compute explicitely basic characteristics of the diffusions with 



generators (67) and (68). For both processes we write, in a concise form, the scale function s as 



and the speed measure m as 



s'{x) 



m{x) 



X 



1 

2|l-a/2 



II 



j,2|a/2 



Obviously the killing measure in both cases is 0. With these explicit formulas at hand we are able to 
show that all boundary points are non-singular (see [5], Ch. II), that is, the process reaches each boundary 
point with positive probability and also starts from this point (at the boundary). We additionally require 
that the speed measure of each boundary point has the value zero, which is consistent with the reflection 
at the boundary. We carry out the appropriate calculation for the point (—1) only. The remaining cases 
can be calculated in the same way. Now, for — 1 < a < z < 1 we have 



m(a, z) s{da) 



2dx 



a (l-a;2)"/2 

dx 

a (1+^ 



da 



(1 - a2)i-"/2 
da 



(l + a)i-"/2 



< 00 
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which shows that the point (—1) is the exit. Similarly, for the same a and z as above we obtain 



{s{z) — s{a)) m{da) 

1 



dx 



(l-x2)l-"/2 

dx 

a (1 + X)W2 



2 da 



(1 - a2)«/2 
da 



:i+a)°/2 



< oo , 



which means that the point (—1) is the entrance. Therefore, the point (—1) is non-singular. Similar 
calculations show that the point 1 is non-singular for the process Xi as well as X2. 
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